Exotic gapless spectrum induced by frustration in quantum antiferromagnets 
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We show strong numerical evidence in favor of an unexpected virtually gapless spectrum, with 
edge states localized at the boundaries, in frustrated spin-1/2 antiferromagnetic ladders with an 
odd number of legs. These features can be accurately reproduced by using a projected BCS wave 
function with a non-trivial pairing, that mixes even and odd reflection symmetries. This approach 
gives the correct classification of the excitations and provides a simple and very appealing picture 
of an unconventional spin-liquid phase stabilized by frustration. 
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In the last few decades there has been an increas- 
ing theoretical and experimental effort to clarify the na- 
ture of the disordered phases stabilized by competing 
interactions in magnetic materials. [1] A very exciting 
scenario appears when such disordered spin liquids are 
gapless with fractional excitations, generalizing the crit- 
ical phase of one-dimensional systems to higher dimen- 
sions. [2J Recent realizations of frustrated antiferromag- 
nets on quasi-two-dimensional lattices, like NiGa2S4 or 
k— (ET) 2 Cu2(CN) 3 , 0, 0] give a promising evidence in 
this direction. In particular, in organic materials, NMR 
shows a power-law behavior of 1/Ti at very low tem- 
peratures without any signal of magnetic order. |5J] In- 
terestingly, also in CS2C11CI4, Neutron Scattering mea- 
surements show the existence of a continuum of excited 
states, compatible with pairs of spinon excitations. [6J] 

From a theoretical point of view, a promising strategy 
is to consider quasi-one-dimensional systems, where sev- 
eral very reliable techniques are available, e.g., bosoniza- 
tion or density-matrix renormalization group (DMRG), 
in order to extract important insight into the more rele- 
vant two-dimensional (2D) case. In this Letter, we take 
this point of view and study the spin-1/2 J\ — Ji Heisen- 
berg model on ladders with an odd number of legs 

n = j 1 

(R,R>) ((R,R')) 

where is the spin operator on site R = (x,y), and 
the sum is restricted to first (Ji), second (J2) nearest 
neighbors. We consider systems with N — L x n sites, 
where n is the odd number of legs and L is the number of 
rungs, with open boundary conditions along the rungs. 

The physical properties of odd-leg ladders are expected 
to be similar to the ones of the 2D case. Indeed, although 
for any finite number of legs n there are no magneti- 
cally ordered phases, in the weakly frustrated regime, 
i.e., J2/J1 <C 1 or J1/J2 <C 1, the excitation spectrum 
is gapless with a well defined spin-wave velocity, anal- 



ogously to the 2D case. Such a system, for small n, 
has been recently considered by different groups with 
contradicting results. @, @, In particular, in Ref. [§| 
it has been argued that, for any number of legs, a rel- 
evant operator that breaks the translational symmetry 
induces spontaneous dimerization in the intermediate re- 
gion J2/J1 ~ 1/2. On the contrary, such a dimerization 
has not yet been detected by numerical calculations. Q 
In order to clarify this issue, we consider the model of 
Eq. |l| with n = 3, by using Lanczos exact diagonal- 
izations and DMRG. Moreover, by means of an accurate 
variational wave function we are able to interpret our 
numerical results. 

Let us start with the excitation spectrum obtained by 
Lanczos diagonalizations with periodic boundary condi- 
tions along the chains. In Fig. [TJ we show the evolution 
of few low-energy excitations in the relevant subspaces 
defined by the spatial symmetries of the Hamiltonian: 
the crystal momentum along the x direction, the reflec- 
tion 1Z X that changes x — > —x, the reflection lZ y across 
the central chain, and the total spin. For both small and 
large J2, as expected, all the excitations odd under lZ y 
(i.e., r y = —1) have a sizable gap. Although the quantum 
numbers of the ground state do not change by varying 
J2/J1, we have evidence in favor of an avoided crossing, 
that is indicated by the sudden drop of the lowest energy 
gap in the ground-state subspace. This suggests the pos- 
sible occurrence of a first-order phase transition, slightly 
shifted to larger values of J2/J1 (i.e., J2/J1 ~ 0.6) with 
respect to the classical model. [1] More interestingly, the 
frustrating interaction induces a dramatic effect on ex- 
citations with odd 1Z y , and for J2/J1 ~ 0.55 the lowest 
excitation is an r y = —1 singlet with zero momentum 
k x (referenced to the ground state). Although Lanc- 
zos diagonalizations are limited to 12 x 3, we have ev- 
idence that, in the thermodynamic limit, this excitation 
becomes gapless in a small region around J2/J1 ~ 0.55. 
It should be emphasized that the quantum numbers of 
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FIG. 1: (Color online) Some of the relevant low-energy gaps as 
a function of J2/J1 from Lanczos data. Triplet with k x = tt, 
r x = 1, and r y = 1 (blue circles), singlet with k x — tt, r x — —1 
and r y — 1 (green triangles), singlet with k x = 0, r x = 1, and 
r a = —1 (red squares), singlet with k x = tt, r x — 1, and 
r y — —1 (empty black triangles), and the lowest excitation 
in the same subspace of the ground state (empty magenta 
circles). All the quantum numbers are referred to the ground 
state and lines are guides to the eye. 



this excitation are not the ones implied by the Affleck- 
Lieb-Schulz-Mattis theorem, i.e., a gapless state with mo- 
mentum k x — tt with respect to the ground state. (loj In 
agreement with that, we have also evidence that both 
a triplet state (r x — 1 and r y = 1) and a singlet state 
(r x = —1 and r y = 1) with k x = tt become gapless 
in the thermodynamic limit for all the values of J2/J1 
considered. These results are not compatible with con- 
ventional dimerization, where the ground state is doubly 
degenerate with a finite gap with respect to all the other 
excitations. 

A further confirmation of a gapless triplet excita- 
tion comes from a systematic finite-size scaling by using 
DMRG with open boundary conditions in both x and 
y directions. The triplet gap vanishes in the thermody- 
namic limit for all the values of J2/J1 considered (see 
Fig. [2]), in agreement with previous DMRG results. Q 
In order to corroborate this conclusion we have stud- 
ied the influence of a third-nearest-neighbor interaction, 
J3. As suggested in Ref. [ill, for J 2 = this stabilizes 
a phase with a sizable gap and dimerization. In this re- 
spect, we have studied the size scaling of the spin gap for 
J2/J1 = 0.55 for different values of J3/J1. As shown in 
Fig. [21 the thermodynamic limit of the spin gap is clearly 
finite for J3/J1 > 0.2, and vanishes as J3 decreases. Since 
the spin gap is characterized by a Kosterlitz-Thouless 
(KT) behavior, the precise value of the critical coupling 



FIG. 2: (Color online) Size scaling for the dimer susceptibility 
X (upper panels) and for the triplet gap A (lower panels) 
evaluated by DMRG. 



J| is difficult to assess numerically. However, our results 
indicate that for J3/J1 < 0.1 the spin gap is exceedingly 
small (< 10 _8 Ji). In addition, extensive calculations on 
five- and seven-leg ladders indicate that the finite-size gap 
at fixed L decreases with the number of legs, thus sug- 
gesting the stabilization of a gapless and homogeneous 
phase for J2/J1 — 0.55 also in 2D. This conclusion is 
also confirmed by a study of the dimer susceptibility cal- 
culated by adding a small perturbation that breaks the 
translational symmetry, O = 5^2 R e l< ^ R Sji ■ Sr+ x , with 
Q = (tt, 0), and by computing the second derivative of 
the ground-state energy with respect to S (see Fig. [5]). A 
clear dimerization is found for J3/J1 = 0.5 but not for 
J3 = 0, as in the latter case the dimer susceptibility does 
not display the sharp \ oc A^ 2 divergence required in a 
symmetry broken phase. [121 ] 

In the following, we will show that a simple variational 
ansatz is able to explain the gross features of the anoma- 
lous low-energy spectrum. Our construction is based on 
a projected BCS (p-BCS) wave function: 



\p-BCS) =V\BCS), 



(2) 



where V is the projector onto the subspace of single oc- 
cupied sites, and \BCS) is the ground state of 



n 



BCS 



= -t 



(R,R')<7 



R,R' 



(3) 

where c' R a (c_r, ct ) creates (destroys) an electron of spin a 
at site R = (x, y), t = 1 is the nearest- neighbor hopping 
amplitude, and Arui = A^/^ are real singlet pairing de- 
termining the symmetry of the BCS order parameter. If 
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FIG. 3: (Color online) Triplet gap Aoo as a function 
of J3/J1 in the thermodynamic limit, from DMRG data. 
The line is a fit according to the KT theory, A 00 = 
Aexp[-b/^/(J 3 - J|)/Ji]. Inset: Detail for J3/J1 < 0.3. 

the pairing involves only nearest neighbors sites along the 
coordinate directions, A x and A v (even under both re- 
flections 1Z X and TZ y ), the BCS eigenstates can be labeled 
by further quantum numbers, i.e., ph = ±1, associated 
to the particle-hole transformation 

PhC^Ph = sign(<j)(-ir+yc R ,- a , (4) 

which commutes with both TLbcs and V, and the pro- 
jected state ([2]) has the same spatial symmetries of the 
spin Hamiltonian. However, the physical states that sur- 
vive after projection are only those with ph = 1 for even 
S z or ph = — 1 for odd S z . Therefore, all physical states 
have a given Ph according to their z component of the 
total spin. In the BCS Hamiltonian, we have the free- 
dom to take either periodic boundary conditions (PBC), 
i.e., c x>y = c x +L,y, or antiperiodic ones (APBC), i.e., 
c x _ y = —c x+ L, y . Of course, in both cases, after projec- 
tion the wave function @ describes a spin state consis- 
tent with PBC. The lowest-energy state is obtained with 
PBC for L = 4m + 2 and with APBC for L = 4m, while 
the lowest excitations correspond to the Gutzwiller pro- 
jection of the BCS ground states with the other choice 
of the boundary conditions, namely APBC (PBC) for 
L = 4m+2 [L = 4m). In this case, the ground state of ^ 
is degenerate because of the presence of four zero-energy 
Bogoliubov modes, that can be identified as spinons, car- 
rying spin 1/2 and momenta k x = ±7r/2. With these 
objects, we can construct four zero-energy states with 
Ph = 1 (that survive after projection): One triplet and 
three singlets. Indeed, two spinons can form i) a triplet 
of momentum k x = 7r even under 1Z X , ii) one singlet with 
k x = ir odd in TZ X , and Hi) one singlet with k x = even 
under 1Z X . A further singlet can be obtained by com- 
bining zero and four spinons, yielding iv) a state with 
momentum k x = n and odd with respect to 1Z X . All 



TABLE I: The low-energy states predicted by the p-BCS 
approach compared with the exact diagonalization data for 
J2/J1 = 0.55 on the 10 x 3 ladder. The variational parame- 
ters are A x = 1.438, A v = -0.987 and A xy = 0.457. 
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these states have lZ y even. The two singlet states ii) and 
iv), belonging to the same symmetry subspace, do not 
represent distinct excitations: In fact, numerical calcula- 
tions show that their overlap increases with the size of 
the system. Analogously, the singlet Hi) is asymptot- 
ically identical to the ground state. Therefore, out of 
the four distinct states, we just obtain two independent 
excitations: A triplet and a singlet, precisely reproduc- 
ing the lowest level of the "tower of states" predicted by 
conformal field theory in the non-frustrated case. [l3| 

Much more interesting is the case when the BCS 
Hamiltonian breaks some symmetry that is instead re- 
stored after projection. For instance, as already empha- 
sized in 2D, [14j we can add a next-nearest neighbor pair- 
ing A xy , with odd reflection symmetry under 1Z X and lZ y . 
In this case, both reflection and particle-hole symmetry 
(|4"|) do not commute with Hbcs since the simultaneous 
presence of A x , A v , and A xy . Instead, after projection 
onto the physical subspace with singly occupied sites, 
these symmetries are restored. Indeed, in this case, the 
BCS Hamiltonian is invariant under 1Z X = Ph 1Z X and 
li y = PhTty Therefore, the ground state of TLbcs has 
a well defined value of 1Z X and lZ y and, after projection, 
since Ph is the identity in the physical Hilbert space, it 
has also well defined values of the true reflection sym- 
metries 1Z X and TZ y . Since Ph is no longer a symmetry 
of Ti-BCSi the eigenstates that were previously forbidden 
in the physical spectrum due to their wrong value of Ph 
are instead now allowed. This approach predicts to a 
pair of new r y = — 1 excitations not present in the usual 
bosonization analysis: They are both singlets, even with 
respect to 1Z X and have different momenta: v) k x = 
and vi) k x = tt. Remarkably, the first of these states has 
precisely the same quantum numbers of the lowest exci- 
tation emerging in the frustrated regime J2/J1 ~ 0.55, 
while the second one, although higher in energy, has a 
sudden drop in the same region, see Fig. [T] 

Further insight on the physics underlying this result 
comes from the structure of the BCS excitations. The 
key feature induced by a non-vanishing A xy coupling is 
the localization of the spinon wave function tpj near the 
edges of the ladder (here j labels the number of the leg): 
ip 2 j = and tp2 j+1 = with r = -(t-2A^ + zA^)/(i + 
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FIG. 4: (Color online) "Phase diagram" of the J1 — J2 three- 
leg ladder found by the variational p-BCS state (see text). 
The Brillouin zone for the 2D lattice is shown in the inset: 
Stars indicate 2D gapless spinon excitations for A xy — [at 
(±7r/2, ±tt/2)], red dots for A xy 0. 



2A xy + iA v ). Only in the A xy -> limit we get r = -1 
and the spinon delocalizes among the chains. Instead, 
if |t| ^ 1, the gapless spinons become edge states. This 
scenario is confirmed by a DMRG study on three- and 
five-leg ladder with odd number of sites per leg, so to 
have a free spinon trapped in the lattice. 

After this detailed exam of the general properties of 
the p-BCS wave function, we show in Fig. [3] the "phase 
diagram" of the three-leg ladder, obtained by variational 
Monte Carlo. The optimal couplings A x , A y and A xy 
have been numerically determined. The most relevant 
feature is the stabilization of a non-vanishing A xy for 
J2/J1 ^ 0.44. The variational estimate of the tran- 
sition point is likely to be somewhat inaccurate, e.g., 
for one chain, spontaneous dimerization is found for 
J2/J1 ^ 0.15 instead of J2/J1 ~ 0.24. Nevertheless we 
believe that the global trends emerging from this phase 
diagram are representative of the real behavior of the 
model. Furthermore, we checked that dimer correlations 
do not display any clear tendency towards ordering. Note 
that, whenever the spinon spectrum is gapped (e.g., when 
PBC are taken along the rungs), the same p-BCS wave 
function describes a dimerized state for finite n. [l5| 

Finally, by considering the set of variational parame- 
ters obtained by Monte Carlo optimization of the ground- 
state wave function, we evaluated the elementary excita- 
tions on the 10 x 3 by a change in the boundary conditions 
along the chain, without a further optimization. In this 
way, beside the ground state, we obtained an explicit 
form of the excited states i), ii), v), and vi) discussed 
previously. The variational energies and the overlaps 
with the exact lowest eigenstates are shown in Table Q] 



for J%l J\ — 0.55. The impressive accuracy of the p-BCS 
wave function for the full set of low-lying states provides 
a clear evidence in favor of the variational picture. 

In conclusion, our results point toward the existence of 
gapless phases in the three-leg spin- 1/2 J\ — J2 Heisen- 
berg antiferromagnet for all values of the next-nearest 
neighbor frustrating interaction. An unconventional gap- 
less phase with an excitation spectrum characterized by 
the presence of low-energy edge states is stabilized upon 
frustration. Although we cannot exclude the possibility 
of an exponentially small dimerization (e.g., with a spin 
gap for three legs of order 10 -8 J\), we believe that the 
new qualitative features that we found faithfully char- 
acterize the physics of the system. This homogeneous 
phase is accurately described by a projected BCS wave 
function representing an algebraic spin liquid. This sce- 
nario is even more plausible in 2D, especially consider- 
ing that the exponentially small tail of the gap implied 
by the KT behavior should disappear, being replaced by 
a more conventional power-law with Aoo ~ (J3 — J3)" , 
and that the same type of wave function remains accu- 
rate when increasing the (odd) number n of legs. The 
strongly frustrated region in 2D is then a gapless state 
with incommensurate spinon excitations, that naturally 
result from the stabilization of a finite A xy , [lj] as shown 



in Fig. [H Such an exotic spin spectrum can be experi- 
mentally detected in 2D frustrated antiferromagnets, like 
Li2VOSiC>4 and VOM0O4 under pressure. 
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